Sandpile behavior in discrete water-wave turbulence. 
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Abstract: I construct a sandpile model for evolution of the energy spectrum of the water waves in 
. ■ finite basins. This model takes into account loss of resonant wave interactions in discrete Fourier space 
• S ■ and restoration of these interactions at larger nonlinearity levels. For weak forcing, the waveaction 
P5 ■ spectrum takes a critical uj~ 10 shape where the nonlinear resonance broadening overcomes the effect 
i of the Fourier grid spacing. The energy cascade in this case takes form of rare weak avalanches on 
"xl" \ the critical slope background. For larger forcing, this regime is replaced by a continuous cascade and 
^ ■ Zakharov-Filonenko uj~ 8 waveaction spectrum. For intermediate forcing levels, both scalings will be 
i/S ■ relevant, u~ 10 at small and u~ s at large frequencies, with a transitional region in between characterised 
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by strong avalanches. 

Keywords: Water waves, wave turbulence, sandpile models, four-wave resonance, energy and wave- 
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^3 ; 1 Introduction 

• i—i . Importance of the basin finiteness for the statistical evolution of the free water surface was recently 
^ \ argued in [1,2] and [3]. Later in the present paper, we will derive an estimate (15) according to which 
c3 ' water surface waves of steepness a are sensitive to the basin size L if their wavelength A is of the order 
or greater than La 4 . This means that waves with steepness a ~ 0.1 and 1 meter wavelength will "feel" 
the boundaries for for lakes or gulfs up to 10 km wide. 

Recent direct numerical simulations of the free water surface in finite basin in presence of gravity 
revealed a bursty character of the nonlinear energy cascade from small to large wavenumbers [2] . Namely, 
cascade strengths were measured as functions of time at two different wavenumbers within the inertial 
range. Intermittent bursts were observed on these graphs, initially arising at the lower wavenumber 
and then propagating to the higher wavenumber. This behavior reminds sandpile avalanches moving 
from low to high wavenumbers and their mechanism can be understood as a cycle: 
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• A cascade arrest due to the wavenumber discreteness leads to accumulation of energy near the 
forcing scale. 



• This leads to widening of the nonlinear resonance. 

• Sufficient resonance widening triggers the cascade thereby draining the turbulence levels and 
returning the system to the beginning of the cycle. 

Note that presence of forcing is essential for this scenario and it should not be expected in freely decaying 
fields, e.g. like in [3]. Below, we present and study a simple model for such a sandpile-like evolution of 
the water wave spectrum. 

2 Differential approximations for waves on an infinite surface. 

Evolution of random weakly nonlinear gravity surface waves in basins of infinite size and depth is well 
described by the Hasselmann kinetic equation [4]. However, for many purposes one can use a simpler 
differential equation model which preserves many properties of the Hasselmann equation [5-7]: 

h — — ^ n 4 jj 26 d 2 \ q\ 
where C\ is a dimensionless constant. This equation preserves energy 

E = ^ J uj 4 nduj (2) 

and the waveaction 

2-7T f 

N = — cu^ndcu. (3) 
Q J 

Even simpler differential approximation was suggested in [7], 

where C<i is a dimensionless constant. This equation conserved both energy and the waveaction, but it 
does not have thermodynamic solutions corresponding to equipartition of these quantities. 

If we insist that having thermodynamics equilibria is important, but ignore conservation of the 
waveaction (which can be done at scales smaller, but not larger, than the forcing scale) then we can 
use another 2nd order differential equation, 

C d ( 2 24 & . S\ . s 
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where C is a dimensionless constant. Similar approach in Navier-Stokes turbulence is called Leith model 
[8,9]. Like in Leith model [9], we can now find a "warm cascade" solution, i.e. the general stationary 
solution of (5) which contains both finite flux and finite temperature components, 



n = — 



I ^-+7- I (6, 



where P and T are (dimensional) constants measuring the energy flux and the temperature respectively. 
For T = we recover the pure cascade Zakharov-Filonenko state, 

n=(P/7C) 1/3 ^ 2 c- 8 , (7) 

and for P = we get the pure thermodynamic distribution, 

n= T -. (8) 



3 Finite-basin effects. 

Let us now consider waves in a square basin with sides of length 2n, so that the wavenumbers take 
values on a discrete lattice, k G Z 2 . The main effect of the finite basin size is in loss of wavenumber 
resonances due to the wavenumber discreteness [1-3]. Indeed, let us consider the 4- wave resonance 
conditions, 

ki + k 2 = k 3 + k 4 , (9) 
cj(k 1 )+cj(k 2 ) = cu(k 3 ) +cu(k 4 ). (10) 

Two different classes of such solutions were found in [2]: collinear quartets (all four wavevectors are 
parallel to each other) and "tridents" (first two wavevectors are anti-parallel and the other two are mirror 
symmetric with respect to the direction of the first two). Parametrisation of the collinear quartets and 
quintets was done in [16]. Note that the collinear quartets are physically unimportant because of the 
zero nonlinear coefficient for such wavevectors, but the next order (5-wave) is nontrivial and yields an 
interesting kinetic equation [16]. The second class, tridents, can be paramertised as follows [2], 

k=(a,0), ki = (-6,0), k 2 = {c,d), k 3 = (c,-d) 

with 

a= (l 2 + m 2 + lm) 2 , b = (I 2 + m 2 - Imf , c = 2lm{l 2 + m 2 ), d = l A - m 4 , 

where / and m are integers. New solutions can be obtained by further rescaling and rotating these 
tridents by rational angles. 

Further, more resonances appear due to the nonlinear resonance broadening even when this broad- 
ening is significantly less than the wavenumber grid spacing [2]. However, the total number of both 
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exact and approximate resonances remain significantly depleted with respect to the continuous case 
and, therefore, they are inefficient for supporting the turbulent cascade unless the nonlinear resonance 
broadening becomes of order of the fc-grid spacing. This allows us to formulate a simplified model for 
wave turbulence in finite basins as explained in the next section. 



4 Wave turbulence in finite basins. 

First, we need to evaluate the 4- wave resonance broadening which is of order of the characteristic 
nonlinear time t^l- Estimate for t^l will be the same if one finds it using (1), (4), (5) or the original 
Hasselmann's kinetic equation [4]. It can also be found from a simple dimensional argument and the 
result is 

r NL ~ s 10 ^ 19 n- 2 . (11) 
This corresponds to the resonance broadening in the /c-space given by 

n NL = 1 J—~g-"J»n\ (12) 

In our model, we will postulate that the wave spectrum will not evolve at uj if the resonance broadening 
unl is less that the fc-grid spacing k, and it will evolve as in the continuous case for knl > k. Such a 
"frozen turbulence" state was first observed in numerical simulations of the capillary waves [11] and it 
was later discussed in [2,12]. One has to be careful, however, not to interpret literally the absence of 
the spectrum evolution at small amplitudes because a small number of exact resonances do survive for 
the gravity waver waves (see the previous section) and further resonances may re-appear at resonance 
broadening which is much less than the /c-grid spacing [2]. However, the number of such resonant 
modes is too small to evolve the spectrum efficiently and in our simple model we just put the Heaviside 
step function H(k nl — k) as a pre-factor to the equation (5) in order to get a model for the spectrum 
evolution in discrete /c-space, 

n = g3/2u4 — [n u, —(um) ) + >y(u)n. (13) 



Here, we have added a function 7(0;) which models forcing at low cu's (e.g. by wind) and dissipation 
at high u;'s (by wavebreaking) and which, in principle, can be a function of n. We will assume that 
in between of the forcing and dissipation scales there exists an inertial range where 7^0. From now 
on, we ignore the dimensionless order-one pre-factor C since our resonance broadening is given by an 
order-of-magnitude estimate. 



5 Behavior predicted by the model. 

Equation (13) with k^l given by (12) will be our master model for the water-wave turbulence spectrum 
in a finite basin. Let us qualitatively consider the consequences of this model. Let us assume that 
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initially there is no waves in the basin and let us start forcing the system at low frequencies. Then, 
there will be no transfer over scales and the spectrum will grow with the growth rate 7 until it reaches the 
critical value where knl = k. After that the nonlinear transfer will get activated and the energy will spill 
into the adjacent range of frequencies. If the forcing is so weak that its characteristic rp = l/max{7(u;)} 
is much longer than the characteristic nonlinear time tnl then the level of turbulence will never greatly 
exceed its critical value and the critical spectrum with will gradually occupy the entire inertial 

range. Condition k nl k gives for the critical spectrum 

n c ~ g 11 / 2 ^- 10 . (14) 

It is useful to re-write this relation in terms of the water surface angle a characterising the wave 
steepness, 

a c ~ (A/L) 1//4 , (15) 

where L is the box size and A is the wavelength. One can interpret this relation as an expression for 
the minimal steepness for which the finite box effect can be ignored. For example, for a box containing 
10 4 wavelengths the finite box effects can be ignored only for a > 0.1. 

If the forcing is stochastic then the subsequent evolution will consist of small avalanches going down 
the critical slope with time intervals At greater than the time ~ tnl needed (according to (13)) for 
the avalanche to travel from the forcing to the dissipation scale. Note that this is a classical condition 
for sandpile models. To be specific, let us consider a type of forcing such that after each interval At 
we add an increment (Aa)e^ where phase is random and uniform in (0, 2n] and a is a small positive 
value, Aa <C y/n c . If at some moment of time n = n c at the forcing scales k G (kp, kp + Ak) then 
with probability 1/2 the spectrum will get greater than critical in the forcing range after time interval 
At and approximately n = n c + Aa^/n c . For At ^> t^tl, such disturbance will have enough time to 
travel /diffuse away before the next spectrum disturbance might appear at k F after another At interval. 
Thus, the evolution of each super-critical disturbance can be treated separately. Because each of such 
disturbances is small, one can use the linearised version of the evolution equation (13), 

with initial condition 

n\ t= o = Aa\/n c (uF) for u G uof + Aa;, and n = otherwise. (17) 
Equation (16) can be re-written as 



^(4^ + ^), (18) 



de d 2 e 

where e = un is the spectral energy density. According to this equation, the disturbance generated by 
forcing will propagate toward high k with speed Ang 19 ^ 2 while getting diffused at an increased rate (due 
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to moving to higher u;'s). The stationary solution of (18) decays as e ~ u~ 3 orn ~ co>~ 4 , i.e. significantly 
slower than n c ~ cj~ 10 . Therefore, for a long enough inertial range the linear approximation will fail at 
some u = u* and the critical slope n ~ uj~ 10 will be replaces by the Zakharov-Filonenko slope n ~ uj~ 8 
for ijj > uj*. The transitional range with u ~ c<j* will be characterised by strong avalanches. 

For stronger forcing, transition to the Zakharov-Filonenko spectrum occurs at lower frequencies or 
even right at the forcing scale if u F > u* (i.e. when a at the forcing scale is steeper than a c at this scale). 
In numerical simulations, efforts are typically made to overcome "frozen turbulence" and generate the 
cascade. At the present level of resolution (up to 512 2 modes) this goal can be achieved with only 
partial success because, according to estimate (15), the condition that turbulence is not frozen can be 
only be marginally reconciled with the condition for the Wave Turbulence theory to work, a < 1. Thus, 
in all existing simulations (e.g. [1-3, 13-15]) turbulence, although not frozen, was still quite sensitive 
to the finite box effects. This state was named mesoscopic turbulence in [3]. In presence of forcing, it 
shows up via strong cascade avalanches coexisting with Zakharov-Filonenko state occupying about a 
decade long wavenumber interval. 

6 Discussion 

In this paper, we presented an evolution model for the spectrum of gravity water waves in finite basins. 
It has the following features: 

• The model is give by a nonlinear second order equation in Fourier space. 

• The model has the cascade Zakharov-Filonenko spectrum and the thermodynamic spectrum 
among its solutions. It also has a general stationary solution where both the flux and the tem- 
perature effects are present. 

• The model takes into account the /c-space discreteness by switching off the nonlinear evolution 
when the spectrum falls below a critical value. The critical spectrum is determined by the condi- 
tion that the nonlinear resonance widening is equal to the /c-grid spacing. 

Based on this model, we established that for very weak forcing the spectrum takes the critical 
slope n ~ uj~ 10 with occasional weak "avalanches" running down this slope. For larger forcing the 
system does not feel discreteness and the spectrum takes the Zakharov-Filonenko form, n ~ uj~ 8 . For 
intermediate levels of forcing, the spectrum may have the —10 exponent at low frequencies and the 
—8 at large frequencies within the inertial range. Such intermediate case is characterised by strong 
avalanches down the mean spectral slope, a feature observed in recent numerical simulations of the free 
water surface [2]. 

It is interesting that steeper than Kolomogorov slopes were also previously obtained for waves with 
narrow-band forcing [17]. The narrow band forcing also leads to quasi-discrete character of the mode 
excitations. To conclude, it is worth mentioning that there are plenty of other examples of dispersive 
waves whose resonant interaction may be affected by the finite size box and where one could expect 
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similar avalanche-like behaviour. Interestingly, irrespective to discreteness, the sandpile analogy have 
also been previously invoked in the wave turbulence context to illustrate sudden readjustments necessary 
to balance the turbulence sources and sinks [18]. 
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